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Abstract. The recent discovery of a diffuse neutrino flux up to PeV energies raises 
the question of which populations of astrophysical sources contribute to this diffuse sig¬ 
nal. One extragalactic candidate source population to produce high-energy neutrinos are 
Blazars. We present results from a likelihood analysis searching for cumulative neu¬ 
trino emission from Blazar populations selected with the 2nd Fermi-LAT AGN catalog 
(2LAC) using an IceCube data set that has been optimized for the detection of individual 
sources. In contrast to previous searches with IceCube, the investigated populations con¬ 
tain up to hundreds of sources, the biggest one being the entire Blazar sample measured 
by the Fermi-LAT. No significant neutrino signal was found from any of these popu¬ 
lations. Some implications of this non-observation for the origin of the observed PeV 
diffuse signal will be discussed. 


1 Introduction 

In 2013, the IceCube collaboration reported on the detection of an astrophysical diffuse neutrino flux 
in the 40TeV - 2 PeV energy range Ul. Since then it has been measured down to about lOTeV in 
energy and the spectrum has been shown to be consistent with a power law of index -2.5 Cl. The 
origin of the observed neutrino signal is so far unidentified. A prominent candidate is the population 
of Blazars, active galactic nuclei (AGN) with their relativistic jets pointing along the line-of-sight. 
Neutrino emission from Blazars has been discussed extensively in hadronic jet models originating 
from p-p or p-y interactions (e.g. ||3) which also produce gamma rays from neutral pion decays. The 
extragalactic non-thermal radiation background at GeV energies and above (EGB) is dominated by the 
emission from Blazars ||4|. An experimental approach is presented here to answer the question if the 
Blazars that dominate the high energy EGB emission also produce a large fraction of the astrophysical 
neutrino flux. The idea is to look for directional clustering in a large sample of muon tracks collected 
by IceCube, around the directions of gamma-ray sources associated with Blazars in the Fermi-2LAC 

cataloglJl- 

This analysis is performed as a maximum likelihood stacking analysis looking for a cumulative 
signal from multiple sources. In contrast to previous stacking searches which looked at 10-30 sources 
0 Q, the populations in this work usually comprise more than 100 sources, the biggest population 
being 862 sources when all Blazars from the 2LAC catalog are combined. Any search for a cumulative 
signal from multiple sources needs to make assumptions about the relative contributions (“weights”) 
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of the individual sources within the population. Here, this weighting scheme is set up in a way to be 
minimally biased with respect to y-v correlation assumptions. Sectionj^covers the motivation for this 
type of search and the technical setup. In Section]^ the sample of muon tracks used in this analysis 
is described. In Section the results of this search are summarized. Conclusions are presented in 
Section|5] 

2 Motivation and Method 

In previous stacking searches for a cumulative neutrino signal from Blazars in IceCube data El 121, 
the measured gamma-ray flux was taken as an estimate for the expected neutrino flux. However, the 
y and v spectra may not be strongly correlated; Other contributions to the gamma flux almost surely 
exist, e.g. inverse Compton photons IQ from high-energy electrons or synchrotron photons from 
muons and protons l9l- Moreover, radiation fields at the sources and in the interstellar medium absorb 
and reprocess the gamma rays originally produced together with the neutrinos. In order to overcome 
this bias as much as possible, two modifications with respect to previous IceCube stacking searches 
have been implemented: 

• There is no sub-selection of sources within the catalog based on their gamma-ray flux 

The populations tested for a neutrino signal are only defined via their spectral classification, the 
largest one being just all Blazars of the 2LAC catalog. Due to the competing production processes 
and different environmental conditions described above, a source with a relatively low gamma-ray 
flux could still contribute significantly to the neutrino signal. 

• Two complementary weighting schemes are used, including a uniform relative weight for all 
sources 

1- ^source ^ 

Note that Fy is the energy flux. Gamma-rays produced in neutral pion decays at multi-TeV 
energies (i.e. by the same process that also produces a neutrino signal) are absorbed quickly by 
IR photons present in the source environment and reprocessed to lower energies. In this process 
they would approximately maintain their bolometric luminosity. 

2 . ^source ~ 1 

Equal weighting is used as the most unbiased weighting scheme, completely ignoring potential 
correlations. Although it is in a strict sense unrealistic, because no population consist of sources 
which all emit equally strong, the resulting limits are quite independent of the unknown corre¬ 
lation between neutrino and gamma-ray flux. They even hold for the case that no correlation 
exists. 

The analysis is performed with an unbinned maximum likelihood Ansatz. The log-likelihood function 
is constructed as 

i 

where N is the total number of neutrino events, 5, and B, are the signal probability distribution func¬ 
tion (PDF) and background PDF evaluations for event i. The background PDF is taken from the data 
distribution in declination and reconstructed energy and is uniform in right ascension. The signal PDF 
is defined as 

nf '"j.tot ■ s j(xi; cri) ■ Sj{Ei) 

^ i ~ iv 9 ^tot ~ ^source ' ^dec. (2) 

V’^'sources ,,, 

Lj 
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where j denotes the individual sources within the population. The term S a i) denotes the point 
spread function (PSF) of source j evaluated at celestial coordinate Xi, depending on the event-specific 
angular error estimate cr,-, while the term Ej{Ei) denotes the energy PDF of source j evaluated for the 
reconstructed energy Ej. The PSF is locally modeled as an analytic 2-d radially symmetric gaussian 
while the energy PDF e is estimated from Monte Carlo simulations. Each individual source within the 
signal PDF is weighted with a declination-dependent detector response weight (Udec. and the relative 
source weight lUsource described above. Equation[T]is minimized with respect to and the likelihood 
ratio with respect to = 0 is calculated for each population and each weighting scheme. It is tested 
if the observed likelihood ratio is compatible with statistical fluctuations of the background events 
by repeating the process for many randomized sky maps. A randomized sky map is generated by 
randomizing the right ascension of the observed events (a more detailed description of this procedure 
is found in ||2l). In case no significant deviations from background fluctuations are observed, upper 
limits on the neutrino flux of the tested population are calculated via the CLj method IfTOl . This is 
performed twice, once for each weighting scheme. 

3 Data and source catalogs 

3.1 Muon Track Data 

A muon-track selection for the years 2009-2011 (IC59-IC86/1) is used, which is similar to the one 
described in QQ The median angular resolution of muon-neutrino tracks with this selection is around 
1° at 1 TeV and 0.5° at 100 TeV. One additional cut is applied which keeps only events that have an 
estimated angular resolution better than 5 degrees. In the southern hemisphere (downgoing events 
in IceCube) this sample is dominated by muons from air showers, while in the upgoing region the 
majority of events come from atmospheric muon-neutrinos. 

3.2 Populations tested 

The populations to be tested are based on the 2nd Eermi AGN catalog (2LAC) Q and are classified 
purely by spectral properties. Only sources that do not suffer from potential confusion are considered 
(CLEAN flag set to True). The populations are listed in table It should be noted that some of the 
populations are overlapping and the p-values correlated (the largest overlap exists between the ESP 
and ESRQ populations, which is about 60 %). 


4 Results 

Eor each population two tests were performed, one for each of the two different source weighting 
schemes. Table summarizes the result^ All test outcomes are compatible with statistical fluctua¬ 
tions of the background. The smallest p-value is 6 % for the "All 2LAC Blazar" sample when tested 
with the "equal weighting" scheme. Eigure shows the resulting neutrino flux upper limits for an 
E^^-^ spectrum for the "All 2LAC-Blazars" population and for the ESRQ population. Notice that the 
flux upper limit using the equal-weighting scheme is represented as a band, instead of a line. It serves 
as an indicator of how robust this flux limit is towards changes of the actual source count distribution 
over the sky, which is unknown. Random realizations are created assuming a source count distribu¬ 
tion for the neutrino emission that has the same shape as the source count distribution in gamma rays 


*The IC79 dataset is the reprocessed one also used in the Blazar stacking search in (TJ. 
^No trial-correction for testing multiple populations is included. 
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Table 1: Definitions of Blazar populations 


Type 

No. of sources 

Motivation 

All 2LAC Blazars 

862 

No bias 

FSRQ' 

310 

BLR radiation ifTTIl 

LSP2 

308 

FSRQ and LSP-BLLAC might be intrinsically similar lfT2l 

ISP 1 HSP2 

301 

HSP objects seem to evolve differently lfT3l 

LSP & BLLAC^ 

62 

Motivated by work in 0 


1 FSRQ/BL-LAC based on optical line equivalent width 

2 Low/Intermediate/High Synchroton Peaked Object (LSP/ISP/HSP), based on position of syn¬ 
chroton peak 


Table 2; Results of the Blazar population tests for both weighting schemes. 


p-values 


^source ^ ^source 1 


All 2LAC Blazars 

36% 

6% 

FSRQs 

34% 

34% 

LSPs 

36% 

28% 

ISP/HSPs 

>50% 

11% 

LSP-BLLACs 

13% 

7% 




(a) (b) 

Figure 1: Neutrino flux upper limits for an spectrum (blue) compared to diffuse bestfit (black 
solid) from la. a) shows this comparison for the "All 2LAC Blazar" sample, b) for the FSRQs. 
Percentages with arrows denote the fraction with respect to the diffuse flux. 
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(parametrization from Cl). For each source in the population a relative source strength is drawn 
from the source count distribution and assigned to the source. The band encompasses the central 
90% of the flux upper limits found from different such realizations. In comparison to the limits the 
best fit diffuse flux above 10 TeV is shown from ja with a spectral index of -2.46. Assuming an 
isotropic flux and 1/1/1 composition of flavors, the maximal contribution from the Blazar populations 
with respect to the diffuse flux can be read off. This contribution to the diffuse neutrino flux is at most 
around 20 %, even if there is a very weak correlation with the observed gamma flux. For FSRQ’s, 
these numbers are lower as the population size is smaller. 


5 Conclusion 

In summary, the Blazars in the 2LAC catalog or any of the tested sub-samples are not responsible for 
the majority of the observed diffuse astrophysical neutrinos. This is in particular interesting since these 
Blazars dominate the high-energy extragalactic gamma-ray sky. The result is also valid if only a weak 
neutrino flux correlation with the observed gamma ray emission is assumed and expected to get worse 
at most by a factor of two for harder spectral indices (fs ^^H-cutoff scenario). Furthermore, one can 
go beyond the 2LAC sources, since extrapolation of the gamma-ray luminosity function shows that 
the 2LAC Blazars (FSRQ’s) already resolve more than 60 % lITSi (70 % lfT6l ') of the EGB contribution 
from all Blazars (FRSQ’s) in the observable universe. Thus, under the condition that the correlation 
between the y and v energy flux holds (i.e. hadronic processes dominate the gamma-ray emission), 
one could further extend the "y-Lum. weighting" limits to apply to the total corresponding population 
by dividing by 0.6 (0.7). 
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